Domain wall Dynamics of the Ising chain in a Transverse field 
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We show that the dynamics of an Ising spin chain in a transverse field conserves the number of 
domains (strings of down spins in an up-spin background) at discrete times. This enables the de- 
termination of the eigenfunctions of the time-evolution operator, and the dynamics of initial states 
with domains. The transverse magnetization is shown to be identically zero in all sectors with a 
fixed number of domains. For an initial state with a single string of down spins, the local magneti- 
zation, the equal-time and double-time spin-spin correlation functions, are calculated analytically as 
functions of time and the initial string size. The domain size distribution function can be expressed 
as a simple integral involving Bessel functions. 



The dynamics of the order parameter, and its cor- 
relation functions are related to experimental observ- 
ables of interacting systems, related to the relaxation 
and transport processes. For systems with known low- 
temperature correlation functions, the dynamics can be 
extremely difficult to determine. For example, the de- 
termination of correlators for hard-core particles moving 
in one dimension has been done using the inverse scat- 
tering method^. The one-dimensional Ising model in 
a transverse field has a complicated dynamics, and real- 
time correlation functions have been shown to have a 
simple structure near a quantum critical point 2]. 

In this paper, we study the zero-temperature dynam- 
ics of a one-dimensional transverse field Ising model. The 
dynamics does not conserve the spin here, and generates 
domain walls, which can propagate and fluctuate in size. 
The evolution of an initial state with one domain of down 
spins (a string of successive down spins in the background 
of up spins), involves creation of more domains and their 
complicated motion. We show that the number of do- 
mains is conserved at discrete times. As a direct con- 
sequence the transverse magnetization is identically zero 
for any state with a fixed number of domains. The calcu- 
lation of the one-point and two-point correlators can be 
carried out using the domain eigenstates. In particular, 
the time evolution of initial states with exactly one string 
of down spins is presented. 

Let us consider the Hamiltonian (on a ring of N sites) 



H — Hz + Hx 



NK 



(1) 



where sf , sf are the components of the spin-1/2 operator 
at site i, K \s, the Ising interaction strength and A is the 
magnetic field strength in x direction. A constant has 
been added for convenience, which we keep along with 
the Ising interaction term in Hz ■ The above Hamiltonian 
can be diagonaHzedlH, and is extensively studied^! (see 
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for a recent review). Let us write the time evolution 
operator as (setting h — 1) 



(2) 



where we have used the interaction representation, 
Hx{t) = exp {iHzt)Hx exp {—itHz). The time-evolved 
operator sf{t) can be expressed in terms of sf , sf^i, sf_i, 
and we have 



2 sin Kt ^ 2 



Kt 



Kt/2 



where the operator di = (sf_^i -f Now, let us 

consider discrete times t„ = rtr, where r = ^-njK (for 
K ~ O.Ole]/, r ~ 10~^^). This simplifies the time evolu- 
tion operator, which becomes \J — exp —itHz cxp —itHi, 



H 



i = A^.f(i-2. 



(3) 



For these discrete times, the evolution operator Uz has 
no role to play, it gives unity as an eigenvalue acting on 
any state with any number of down spins. The operator 
Hi can only act at sites with one up and one down spin 
on either side, in which case it will flip the spin. The 
polarized state, either all up or all down spins has no 
dynamics, as it is an eigenstate of Hi with eigenvalue 
zero. Hi operates on a state with a string of m down 
spins to generate states with either m-\-l or m — 1 down- 
spin strings. In time, the size of the string can change 
(from TO = 1, a lone down spin, to to = iV — 1, a lone 
up spin) and the center of mass of the spring can move 
in either direction. Similarly the operator acting on a 
state with many strings of down spins, causes them to 
move and fluctuate in size, but the number of strings 
will remain the same. 

Let us consider the subspace of one-string states, 
\l,m > where to = l,iV — 1 denotes the length of the 
down-spin string starting from the site I ~ l,N. By ar- 
ranging the basis states in the order |1,1 > ..|l,Af — 1 > 
,..|A^, 1 > ..\N, N — 1 >, the Hamiltonian has a struc- 
ture oi N X N tri-diagonal matrix with periodic bound- 
ary conditions. The diagonal positions are occupied by 
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a.N— IxN—l tri-diagonal matrix with open bound- 
ary conditions, and upper (lower) diagonal positions are 
occupied by a matrix with nonzero elements only on the 
lower (upper) diagonal positions. All the nonzero ma- 
trix elements are equal to A/2. The cigcnstates of the 
Hamiltonian Hi can be easily constructed as 

\k,q) = AJ^^"'"'-^^'^"'^ sm{qm)\l,m), (4) 

Lrn 

where A is a normalization constant, and the two mo- 
menta q = ttI/N, I = 1, N — 1 and p = nn/N, n = 1,N 
label the one-string eigenstates. The phase corresponds 
to the uniform motion of the center of mass of a string of 
m down spins. The eigenvalue for this state is e{k,p) = 
2Acosfccos(7. 

The eigenf unctions are the same as the two-magnon 
eigenfunctions of the XY model, or the two-fermion 
tight-binding eigenfunctions, labeled by two momenta 

< qi,q2 < 27r, Vij = Ae^P^^~^^^ sin q{j — i), using 

1 ~ l,j = I + m in the above, and the relative mo- 
mentum q ^ qi — 52/2, and the average momentum 
p = qi + q2/2. In the two-magnon eigenfunction the la- 
bels i,j stand for the two down spins, whereas here they 
stand for the first sites of the down-spin and up-spin do- 
mains in the domain eigenfunctions. Also, the range of 
p here is < p < 27r, whereas for the domain eigenfunc- 
tions, we have < fc < tt. So the two ranges < p < tt 
and TTp < 2'K of the two-fermion eigenfunctions map onto 
one range < fc < tt of the domain eigenfunctions, and 
each value of p corresponds to two eigenstates (in terms of 
the fermions, the two states are 4iT42i 1^^' 4ii42T 

|0)). 

Similarly an n-string eigenstate maps onto a 2n-magnon 
eigenstate of XY model (or a 2n-fermion tight-binding 
problem, with n fermions with spin down, and spin con- 
figuration being a Neel configuration). It is easy to check 
that for any one-string eigenstate, the expectation value 
< Si{t) >= 0, for all times. And similarly for eigenstates 
in other sectors. However, the dynamics of an initial state 
with a given spin configuration is quite complex. As we 
shall see below the transverse magnetization is still zero, 
but, however, (sf) can be nonzero. 

Let us consider at i = a state with a string of tq 
down spins, starting from a site io and ending at the site 
Jo = *o + ''o. At later discrete times, the state can be 
written as 

I0W> -E-^'^wi^'j')' (5) 

where the wave function can be written in terms of the 
single-string eigenfunctions. Let us write the two-fermion 
wave function as 

=E^""^'''^*^'o,.oV'*, (6) 

k,p 

which becomes (after taking N — s- 00) 

<pij = e 2 ( +^ " ( Jj_ju (a;) Jj_j(j — Ji-j„Jj-io) (7) 
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FIG. 1: The local magnetization as a function of the position 
r = i — io for an initial state with ro — 71 down spins, for 
various values of time x = 2At. 

where the argument of all the Bessel functions above is 
X = 2At. And = (pij for j > i, whereas 0^ = 4>i^N+j 
for i > j (here the length of the down-spin string is 
N — i -\- j). Though the fermion wave function is the 
same for both cases, but for the domains it is quite dif- 
ferent, j < i configurations occur with negligibly small 
probabilities for a large N. The expectation values of 
operators at time t reduce to sums over Bessel func- 
tions. Let us first consider the transverse magnetization, 
< s'f{t) >=< sf + s- > /2. The raising and lowering 
operators connect basis states that differ by a unit string 
length, and we have 

< S~ >= "^^m-^^-l-l,™ + (f>*+l,r,i'f>^,r,i = 0, (8) 

and similarly the raising operator has zero expectation 
value. Both the components of transverse magnetization 
vanish identically here. For the case of states with more 
number of strings, the argument is as follows. As the 
wave functions are 2n-fermion wave functions, a nonzero 
expectation value for < s~ > would imply an off-diagonal 
long-ranged order; which is prohibited by an argument 
due to Yang viz. the off-diagonal elements of a one- 
particle reduced density matrix cannot exhibit long range 
order. For a nonvanishing transverse magnetization, one 
has to construct superpositions of states drawn from 
different domain-number sectors. These states will be 
more meaningful for the low-temperature behaviour for 
A> K. 

Let us turn to the magnetization at a given site, and 
since sf is diagonal in this basis, we can express its ex- 
pectation value as, 

i-l i N i N 

l+l i+l l+l i+1 

Let us denote i = io + r, jo = *o + ^Oi a-nd choose io = 
N/2 for convenience. Since r = —a,r = ro — 1 -I- a are 
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FIG. 2: The diagonal cqual-timo correlation function is plot- 
ted as a function of the separation, for a few values of {r,x). 
The initial domain size is ro = 30 in all the cases. 

equidistant locations from the string on either side of the 
string, those two sites will have the same magnetization 
for all times. Thus, it suffices to study r > 0. We have, 

{S'r) = l{hf2-f) (10) 

where /i = g{r), /2 = g{r - ro), for r > ro, /2 = -^(ro - 
r— l),forO < r < ro, and the function 5 (r) = YT-r "^ni^)- 
It can be expressed as 

g{r) = 1-2 Jr{xi)Jr+i{xi)dxi. (11) 
Jo 

The function /(r, ro) = 2^^^^ J„J„_ro can be written 
as 

/(r, ro) = / {JrJr+l-ro + Jr+lJr-ro)dxi 

Jo 

X 

= — {JrJr+l-ro ~ Jr+lJr-roi^)) (12) 

ro 

For a fixed value of x as r — » oo, 5 ^ 1, / ^ 0, which 
implies < >^ 1/2. This is because the domain wall 
will take x ~ ro +r time to reach the far away point. For 
fixed r, as x — > 00, the two functions have the behaviour 
(J ^ 0, / — > 2/7rro, for ro odd, and it is zero when ro is 
even. This implies all the spins within a size of x from 
the initial position have an expectation value 2/7r^rQ 
for ro odd, whereas for even initial string sizes, the lo- 
cal magnetization fluctuates around a zero expectation 
value. The local magnetization is plotted as a function 
of the site index r, for various values of time in Fig.l, for 
an intial state with a string of 71 down spins. 

The even-odd effect can be understood as follows. First 
we note that for ro odd (even), at later times odd (even) 
sizes are more probable than even (odd). Later we shall 
compute the size distribution function. The down-spin 
density is maximum at the center of the domain, and as 



we move away from the center the spin density drops. 
Most of the contribution for a negative magnetization 
should be coming from the probability of a site being a 
center of a domain. The center of mass of even-size do- 
mains cannot be a lattice site. Hence, for even ro, the 
magnetization fluctuates around a zero mean. However, 
this even-odd discrepancy disappears for larger initial 
string sizes, viz. for ro = 9 the mean drops to —0.002 al- 
ready. This is because as ro becomes large, smaller-sized 
domain are less probable, as we shall see later. And large 
odd-size domain move slower, thus depleting the proba- 
bility for a site to be a center of a domain. 

Now we turn to the calculation of equal-time correla- 
tion functions. The off-diagonal correlation function is 
straightforwardly written as (sfsj) = The 
diagonal correlation function can be written as 

j N i j 

(4(t)s|(t)) = i - ^(EE+EDi-^'^i'- (13) 

i+l j+1 1 i+1 

Let us choose, i = io + r,j = i + s,jo = io + ro, and 
io = N/2 as before. Here there are three regimes, de- 
pending whether both i, j are within the initial string, or 
both outside, or one of them is outside. The correlation 
function can be reexpressed as 

(s^rSr+s) ^ \{FlF2 - F') (14) 

where Fi = g{r + s) — g{r) and F = f{r, ro) — f{r + s, ro). 
The function F2 has different dependencies in the three 
regimes, as, F2 = g{r — ro + s)—g{r — ro), for r > ro,-F2 = 
g{ro - r - 1) - ,g(ro - r - s - 1), for r -|- s < ro, = 
g{ro - r - 1) + g{r + s - ro), for r < ro,r -|- s > ro- 
For a fixed value of the separation s ^ 00, these 
functions reduce as Fi /i, -F2 f2,F ^ f. And the 
two-point correlation function decouples, and becomes 
(s^Sr+s) ^{Sr) (for a fixed time, the far away spin is 
up, as the domain wall has not had time to propagate 
till that far). In Fig. 2 the correlation fimction has been 
plotted as a function of the separation, s, for a few values 
of r and time, for initial states with one domain of ro = 30 
down spins. 

The double-time correlation functions are not very 
easy to evaluate. Let us define Cij{t2,ti) = 
{(j)o\sf{t2)sl{ti)\(l)o). Setting ti = 0,i2 = t, we have 

a,(t,0) = «(i))(5^^(0)). (15) 

The operator s| can be taken out, as |^o > is its eigen- 
state, and the rest is just (s|(t)). But for an arbitrary 
ti ^ 0, the correlation function can be written as 

Cij{t2,ti) = ^ \(t>lm{xi)\^l'm'{X2)\{Si)l'm'{Sj)lm 

l,m^l'm' 

(J;'_;(X2 — Xi)Jm'~m — Jv -rnJ-m' -l) ■ (16) 

In the above, the argument of all the Bessel functions is 
the difference of time. The double-time correlation func- 
tion is substantially more complicated, and a detailed 
analysis will be presented elsewhere. 
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FIG. 3: The size distribution function is plotted as function 
of the domain size r, for a few values of time (the initial size 
is ro = 10 in all cases). 

We now turn to the domain-size distribution. The 
probability for a size r to occur will depend on the initial 
domain size, and the time. We can write the probabil- 
ity as p{r,ro,x) — J2\4'ij\'^Sj^i+r, which becomes, after 
replacing the delta function by a Fourier integral, 

P--/ e''^'^[cos{<l>ro)h\0,x)~h^{ro,x)]. (17) 

We have defined another auxiliary function involving a 
sum over Bessel functions, 

h{y, x)^\Y, Jk{x)Jk+y{x)e^''^\ = Jy{2x\ sin ||), (18) 

— oo 

where we have used the summation rule for the Bessel 
functions in the last step|2j. After a few more manipula- 
tions, we get 

2 dt 

p^- ^==[Jr+roi2xt) - i-iy^Jr-roi^xt)]'- 

TT Jo Vl - t2 

(19) 



In Fig. 3 we show the distribution function plotted against 
the domain size for a few values of x, for an initial state 
with one domain of down spins of size tq = 10. The 
probability is maximum for r — tq, and as x becomes 
very large all sizes are equally probable modulo an even- 
odd effect. The above distribution becomes simpler for 
the case of = 1 as 

8^2 /.7r/2 

p{r,l,x)^^ Jl(2xcos4)). (20) 

T^X^ Jq 

This distribution function can change substantially for 
initial states with more than one string of down spins. 
For a finite density n of strings, the domain size r is 
limited by the average domain spacing in the initial state. 
The probability is more for domain sizes r < rc 1/n. 

In conclusion, we have explored the domain-wall dy- 
namics for the Ising model in a transverse field in one di- 
mension. The eigenstates of the time evolution operator 
are simple domain eigenstates for discrete times (multi- 
ples of 2TT/K). The transverse magnetization is identi- 
cally zero in all sectors with different number of domains. 
Explicit calculations of the local magnetization and the 
spin-spin correlation functions are presented for initial 
states with one domain of down spins. Though the above 
analysis is valid for large A also, these domain states do 
not belong to the low-energy spectrum here. One can 
do a similar analysis with domain states using the eigen- 
states of sf operators, which are more natural states for 
A > K. Finally we would like to comment on the do- 
mains in two dimensions, where the situation is quite 
different. Firstly, the domains can have various shapes. 
Secondly, though at discrete times, the time evolution 
can be simplified analogously, however, the number of 
domains is not conserved. Another major difference is 
that the domain size is limited, and they do not move 
easily. For a domain of arbitrary shape at some time, 
the domain size at later times is limited by the size of 
the smallest rectangle into which the initial shape can be 
enclosed. 
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